Self-consistent analysis of electron-phonon coupling parameters of graphene 
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We present a self-consistent analysis of the photoemission spectral function A(k,cu) of graphene monolayers 
grown epitaxially on SiC(OOOl). New information derived from spectral intensity anomalies (in addition to 
linewidths and peak positions) confirms that sizeable kinks in the electronic dispersion at the Dirac energy Eo 
and near the Fermi level Ep arise from many-body interactions, not single-particle effects such as substrate 
bonding or extra bands. The relative electron-phonon scattering rate from phonons at different energy scales 
evolves with doping. The electron-phonon coupling strength is extracted and found to be much larger (~ 3.5 — 5 
times) than predicted. 

PACS numbers: 73.21.-b,71.38.-k,73.22.-f, 79.60.-i 



Many-body interactions in epitaxial graphene are interest- 
ing since graphene has been proposed for numerous device ap- 
plications including high power electronics and novel device 
schemes fEIH]. Full characterization of the electron scatter- 
ing lifetime under a variety of conditions is therefore central 
to understanding the properties of such devices. Moreover, 
graphene is an excellent system for investigating theories of 
many-body interactions in two dimensions. First, it is straight- 
forward to prepare high quality films, with intrinsic photoe- 
mission linewidths among the sharpest available for any ma- 
terial UIU. Second, a large change in the carrier density of 
the order of ±10 13 e~/cm 2 can be achieved through applied 
voltage in a sated device |5|,|g], or equivalently through chem- 
ical doping I2D, suggesting a tunability of many -body effects 
for novel devices. 

We have previously shown that the electron scattering life- 
time of holes in n-doped graphene is dominated by a combi- 
nation of ordinary Fermi liquid electron-hole pair (e-h) excita- 
tions, electron-phonon (e-ph) coupling, and electron-plasmon 
(e-pl) coupling fl]. Theoretical work has qualitatively con- 
firmed the e-h and e-pl interpretation. Both Hwang et al. 13] 
and Polini et al. |@1 have shown that the electron-plasmon in- 
teraction leads to an enhancement of the scattering rate, and 
a "7T-band mismatch", where the lower and upper tt bands 
are displaced from each other when projected through their 
crossing point at the Dirac energy En. These effects are 
similar to those seen in experiment, but are not predicted 
in single-particle theories, whereas they naturally arise from 
many-body interactions. Similarly, the dispersion of the bands 
near Ep is heavily modified, a fact which is ascribed conven- 
tionally to e-ph coupling since the effects occur ona~ 200 
meV energy scale, corresponding to the phonon bandwidth of 
graphene | Tol 11, 121. 

Although the evidence for e-ph and e-pl coupling is strong, 
the existence of the e-pl coupling and the quantitative estimate 
of the e-ph coupling strength A derived from the data remain 
controversial. First, as an alternate model to the e-pl cou- 
pling, the bonding of graphene to the substrate induces a gap 



at En 1 1 3l I14I1 . Alternative causes for a dispersion anomaly 
are the presence of defects or quantum size effects for islanded 
films iH El- 

Second, a simple estimate of the coupling st rength fl7ll gav e 
A~ 5 times stronger than detailed calculati onslldll ll[l8l[l9ll . 
This suggests a stronger role than predicted ! 1 8ll20fl for 7r-band 
e-ph coupling in superconductivity of graphite intercalation 
compounds (GICs) like CaC@. But the linear-band estimation 
method 11711 was shown to overstate the coupling ||2lj], and fur- 
thermore finite energy and momentum resolution was cited to 
explain at least part of the discrepancy |12J] . 

Here we provide a self-consistent analysis of the experi- 
mental spectral function which makes no a priori assump- 
tions about the bare bands, but assumes only non-violation 
of causality, approximate particle-hole symmetry, and a weak 
momentum- (k-) dependence of the self energy. Self- 
consistency between model and data is demonstrated by a 
comparison of the scattering rates, dispersion energies, and 
the absolute spectral intensity - the latter having been not usu- 
ally considered. From our analysis of the low doping regime 
(n < 6x 10 13 e~/cm 2 ), we find that: (i) key spectral features 
arise from many-body interactions, (ii) previous estimates of 
large e-ph constant are confirmed, and (iii) the relative e-ph 
coupling strength to multiple phonon modes is shown to be 
strongly doping-dependent. These results are a challenge to 
the present understanding of the carrier lifetimes in graphene. 

The photoemission single particle spectral function is 



A(k,w) 



|ImS(k,o;)| 



(w - cj b (k) - Re£(k, u)Y + (ImS(k, oj)) 2 



(1) 

where u> is the quasiparticle energy, a>b(k) is the bare (un- 
renormalized) band, and S(k, u>) is the complex quasiparticle 
self-energy, whose real and imaginary components are related 
by Hilbert trasformation to satisfy causality. We seek to deter- 
mine E(k, u>) with no knowledge of the bare band Wb(k) or 
the u> > spectral function, and subject to uncertainty due to 
experimental broadening. Once Re£(k, uj) is determined, the 
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FIG. 1 : (a) Experimental and (b) Simulated spectral functions A{k.,u>) for two dopings n = 2.5,3.9 x 10 13 e~ /cm 2 . The dashed lines are E-p 
and the K point of the Brillouin zone, while the solid lines are the self-consistent bare bands, (c) Real and Imaginary parts of the self-energy 
(with 0.02 eV offset), (d) the difference between renormalized and bare bands, and (e) MDC amplitudes for experiment (red, green) and model 
(black) as a function of doping. The lowermost spectra correspond to as-grown graphene (n = 1.2xl0 13 e~/cm 2 ), with doping increasing 
linearly for each subsequent spectrum, up to n = 5.6x 10 13 e~/cm 2 . The heavy lines in (c-e) mark the approximate Dirac energy En- 



e-ph coupling constant A is given by 

A = -dRe£(k,w)/<9fc| w=0 - 



(2) 



In the usual analysis|22j, I23H . each momentum distribu 



tion curve (MDC) A(k,w=const) is fitted to a Lorentzian 
and A(k, ui) is thereby parametrized into three functions: 
the renormalized band dispersion w(k), the Lorentzian width 
W(u>), and A(oj), the amplitude along tcj(k). Provided the 
bare band is linear, Im£(k, ui) is simply proportional to W(ui) 
and the self-energy is easily found. However, if Wb(k) is non- 
linear the problem is more difficult because of the non-trivial 
relationship between W(oj) and ImS(k, uj). 

To solve the general problem, we apply an optimization ap- 
proach similar to Kordyuk et al. 12411 : From initial guesses 
for Wb(k) and S(k, u>), we calculate the simulated spectral 
function A(h, ui) and by MDC analysis parametrize it into 
three functions W(uj), A(u>), £(k). The self energy and bare 
band are iteratively refined in order to minimize the difference 
cD(k)— w(k) subject to the constraint W(lu)=W(oj). Agree- 
ment between simulated and experimental A(u>) is left as a 
final check of the model. To reduce the number of free pa- 
rameters, we use a quadratic bare band. 

The experiments were performed on in situ grown samples 
at beamline 7.0 of the Advanced Light Source. Samples were 
prepared by epitaxial growth on SiC(0001) and doped with K 
atoms as described elsewhere The instrumental broad- 

enings were 25 meV and 0.01 A -1 , and the photon energy 
was 94 eV. The sample temperature was ~20 K. 

Fig. l(a,b) shows two typical experimental spectral func- 
tions A(]s.,uj) for different ?i-dopings of graphene, in com- 
parison to simulated spectral functions A(k, uj). The latter 
were calculated using the optimized bare bands and S(k, ui) 
shown in Fig. l(b,c). We assumed particle-hole symmetry, 



i.e. that ImS(k, ui) is a strictly even function with respect to 
u) = 0, implying RcS(k, 0)— ► 0. The intensity was scaled 
by a function linear with oj to account for the non-uniform 
sensitivity of the electron detector. Comparison of the quan- 
tity (w(k)-Wb(k)) derived from MDC fits of A(k,uj) and 
A(k, uj) (Fig. 1(d)) shows excellent agreement, demonstrating 
self-consistency of our derived S(k, uS), and confirming that 
deviations in the band dispersion indeed are fully described 
by quasiparticle scattering, and not details of the initial-state 
band structure. In particular, the energy gap at Ed lfl3i[l4ll . if 
any, must be much smaller than the band renormalization by 
many-body interactions. 

We stress that in our analysis the simulated and experimen- 
tal spectral functions are treated on an equal footing with re- 
spect to temperature smearing at E-p and experimental broad- 
ening, both of which are included in our simulated A(k, w). 
Therefore these uncertainties are in principle deconvolved out 
of our derived self-energies. Furthermore, we find that the cal- 
culated and experimental A(u>) in excellent agreement (Fig. 
le), despite not being included in our optimization. This is a 
stringent test of the analysis, because as we now show, A(ui) 
is highly sensitive to many -body interactions, perhaps more so 
than the MDC width and peak positions, a fact we can exploit 
to determine new information about e-ph coupling. 

Fig. 2(a) shows simulated spectral functions for a band in- 
teracting with a single Einstein optical mode at u>o=200 meV 
using realistic experimental parameters, including a power- 
law scattering rate to simulate electron-electron interactions. 
The MDC amplitude function A(u>) shows a relatively flat top 
down to wo- Suppose we add additional weak scatterers at 
lower energy scales (Fig. 2(b,c)); how is the MDC analysis 
affected? Casual inspection of the A(k, uj) images shows that 
scarcely any effect on the apparent dispersion or scattering 
rates can be observed with realistic broadening and limited 
statistics. However, the MDC amplitudes A{lo) are strongly 
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FIG. 3: Experimental A(u) for the same samples as in Fig. 1. The 
onset position of the phonon kink at ~ 200 eV is shown by the 
dashed line. 



FIG. 2: Simulated spectral function and corresponding self-energy 
and fitted MDC amplitude A(uj) for (a) a single Einstein phonon 
mode (b) the same, plus an additional Einstein mode at 100 eV of 
20% relative strength, (c) a single mode at u>o = 200 meV with 
additional modes decaying towards to = 0. 



affected, with an easily visible shift of spectral weight from 
ujo towards Ep. This shows that A(lo) can access information 
about weak or even sub-resolution features in the self-energy. 

We believe that the evolution of such 'hidden' features with 
doping are necessary to explain the evolution of A{lo) near 
Ep. Fig. 3 shows the measured A{lo) as a function of doping. 
While all spectra show a break in slope around 200 meV, there 
is clearly an evolution from a merely prominent shoulder be- 
ginning around 200 meV to a flat top below Ep with increased 
doping. The simulations in Fig. 2 suggest that this shoul- 
der is the dominant mode at all dopings, but its strength in- 
creases relative to the lower energy modes with doping. Such 
a relative evolution of the e-ph coupling with different phonon 
modes has not been predicted by theory. 

We close by discussing how the self-energy and experimen- 
tally derived e-ph coupling constant A compares to recent cal- 
culations. The globally optimized self-energies in Fig. 1 do 
not fit perfectly well near Ep, since the quadratic bare band 
does not fit the data perfectly over all energies. To achieve a 
more accurate S(k, lo) near Ep, we have analyzed the data in 
a smaller energy region, shown in Fig. 4(a) for a sample with 
Ep-Ep, ~ .65 eV. 

The simulated A(k, lo) and the the derived self-energy 
E(k, lo) and bare band t^b(k) are shown in Fig. 4(b, d). Good 
agreement between measured and simulated Lo(k), W(lo), and 
A(lo) are demonstrated in Fig. 4(e-g). For comparison, in 
Fig. 4(c) we show the spectral function calculated from the 
predicted e-ph S(k, lo) from Ref. Il 111 . Both the experimental 
and the predicted A(k, lo) were modeled with a .028 eV offset 
to ImS(k, lo) to take care of the background defect scattering 
rate. Although in excellent qualitative agreement, the theoret- 



ical self-energies are much smaller than those measured (see 
comparison in Fig. 4(d)), and hence the kink, scattering rates, 
and the abrupt increase in the spectral amplitude A(lo) are not 
well-reproduced (Fig. 4(e-g)). This conclusion is not affected 
by experimental broadening as proposed in Hill because the 
same discrepancy appears in both the optimized self-energies 
(from which the broadening is essentially deconvolved, see 
Fig. 4(d)) as well as from the apparent self-energy derived 
from both experiment and broadened theory (Fig. 4(e-f)). Fur- 
thermore the disagreement in A(lo) is practically unaffected 
by our experimental broadening. 

The discrepancy in A ranges from a factor of 3.5 to 5 times 
larger than predictions for e-ph coupling 

EElEl, the lower 
estimate found by scaling the predicted real and imaginary 
parts of S(k, lo) to match the experimental functions, and the 
upper estimate by applying Eq.|2]to the two Re£(k, lo) func- 
tions in Fig. 4(d) (25J. If some of the scattering is attributed to 



electron-electron coupling near Ep (which rises as 



2 )H 



then this would act to reduce the derived e-ph coupling con- 
stant, however, the e-e contribution to ImS(k, lo) is predicted 
to be very weak (around .004 eV at lo = 200 meV) lEr3l . and 
it cannot explain the sharp kink at the phonon energy scale 



Why is the apparent e-ph coupling so large? First, it could 
be due to coupling to some other mode (e.g. magnon) at simi- 
lar energy scale, although such a mode has not been observed 
or even predicted. It cannot be due to coupling to substrate 
phonons, because the SiC phonon modes are at lower energies 
than the observed kink. Furthermore there is a carbon-rich, 
graphene-like buffer layer between graphene and SiC 12711 . 
whose coupling to the graphene electrons would also have to 
be anomalously high to explain the results. 

The presence of our dopant K atoms cannot explain the 
large coupling: first, because of the apparent strong coupling 
even in our as-grown films, and second, because K vibrations 
are at too low an energy to explain the kink (although they 
might contribute in principle to the lower-energy modes sug- 
gested in our self-energy in Fig. 4(d)). Also, the very minute 
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FIG. 4: Spectral function of graphene for n = 3.4 x 10 13 e~/cm 2 (a) experimental A(k,u), (b) optimized model A(k, uj), (c) model A(k, ui) 
using self energy by Park et al. (d-g) the self-energies, dispersions, MDC widths, and A{uj) are compared for the results of the MDC analysis 
of A(k, uj) (markers) and the models, with (solid, dashed) lines corresponding to (b, c). 



amount of K atoms (~ .06 per graphene unit cell at our highest 
doping) makes its presence unlikely to cause the effect we see. 
Lastly, there is a small rippling of graphene on SiC 12811 . Al- 
though a local curvature can enhance the e-ph coupling l29ll . 
the curvature in our films is too small to have a significant 
effect. 

Our findings lead to three possibilities, first, that A for 
graphene is anomalously strong, second, that there is anoma- 
lously strong scattering from something other than phonons 
(such as defects, or e-e scattering) on a similar energy scale, 
or third, there is a mutual interplay of interactions that leads 
to an enhanced scattering rate overall. Noting that the peak 
near Ep, in ImE due to plasmon scattering is broadened and 
shifted towards Ep in our experiment relative to calculations 
JUI^], such an interplay is conceivable. 
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